Lower bounds on the absorption probability of beam splitters 
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We derive a lower limit to the amount of absorptive loss present in passive linear optical devices 
such as a beam splitter. We choose a particularly simple beam splitter geometry, a single planar slab 
surrounded by vacuum, which already reveals the important features of the theory. It is shown that, 
using general causality requirements and statistical arguments, the lower bound depends on the 
frequency of the incident light and the transverse resonance frequency of a suitably chosen single- 
resonance model only. For symmetric beam splitters and reasonable assumptions on the resonance 



frequency ut, the lower absorption bound is p n 



PACS numbers: 42.50.Ct, 42.79.Fm, 42.50.Nn 



I. INTRODUCTION 

Passive linear optical elements such as beam split- 
ters are indispensible tools in photonic interferometry. 
Most recently, they became even more important with 
the emergence of ideas to use beam splitters as central 
units in all-optical quantum information processing 0. 
The action of a beam splitter can be understood by as- 
suming that light impinges on a single slab of material 
with a refractive index r\. The differences in refractive 
indices at the faces of the beam splitter cause light to be 
partially reflected from the surface, while other parts are 
transmitted. 

In a mathematical description of it, the beam split- 
ter is mostly assumed to be lossless, at least in the fre- 
quency window associated with the band width of the 
incident light. In this case, the action of a beam splitter 
on two amplitude operators of photons impinging on it 
is given by an element of the unitary group SU(2) which 
reflects energy conservation by splitting into transmitted 
and reflected light Q- If two identical photons in well- 
defined single-photon Fock states impinge on both sides 
of a lossless beam splitter, one would for instance observe 
the well-known Hong-Ou-Mandel quantum interference 
effect 3] with perfect visibility. 

Realistically, however, perfectly lossless beam splitter 
do not exist and cannot exist. Causality and the re- 
sulting Kramers-Kronig relations for the dielectric per- 
mittivity imply that a (real part of the) refractive index 
7](lu) different from unity is always accompanied by an 
imaginary part k(uj) responsible for describing absorp- 
tion and dissipation £| . An additional feature of the so- 
called superconvergence rule p| which follows from the 
Kramers-Kronig relations is that the real part of the re- 
fractive index is unity on average over all frequencies, i.e. 

On the other hand, causality itself does not imme- 
diately imply that the absorption coefficient has to be 



'Electronic address: |s . scheeKBimperial . ac . uk| 



iQ- b {u)/uj T y 



strictly larger than zero for all non-zero frequencies. In 
fact, the Kramers Kronig relations in principle allow for 
a situation in which there is only one 5 function peak 
in the absorption spectrum at some particular frequency 
and vanishing absorption elsewhere. Physically, this im- 
plies an infinitely sharp resonance and thus an infinitely 
long lifetime of the corresponding transition which is not 
realistic cither. Here some statistical arguments come 
to the rescue. For macroscopically large systems a sta- 
tistical description in terms of response functions has to 
be used which leads to a particular form of the absorp- 
tion spectrum. In the linear-response approximation, 
this prevents the absorption function to become zero at 
any point except the origin of the frequency spectrum. 

A causal response of a macroscopic physical quantity 
(such as the dielectric polarization) to an external exci- 
tation (such as an electric field) results in relations be- 
tween the correlation function for the responding quan- 
tities and the imaginary part of the response function. 
Such relations are known as fluctuation-dissipation the- 
orems |3- They have wide applications in all branches 
of statistical physics. Of particular interest to us is that 
fluctuation-dissipation theorems have helped to establish 
limits on physical processes such as frequency stabiliza- 
tion of lasers with cavities Q • Generalizations of these 
theorems have been used to provide bounds on the per- 
formance of solid-state two-qubit gates @ . 

A realistic description of a beam splitter that includes 
absorption and dispersion has been formulated in [|J and 
makes use of the quantum-optical input-output relations 
derived in [To||. It is shown in Q that the action of a 
lossy beam splitter on the amplitude operators of the 
incoming light (plus the amplitude operators of the de- 
vice excitation) corresponds to an element of the group 
SU(4). From there the SU(2)-action of a lossless beam 
splitter follows as a limiting case in which photons and 
device would effectively decouple which, by the above 
arguments, means that there is no beam splitter action 
present at all. 

The consequence of this argument is that any beam 
splitter, no matter how well it is fabricated, must show 
some non-vanishing absorption. In this article we will 
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show how to obtain lower bounds on the amount of ab- 
sorption. We will use a simple model of a single pla- 
nar slab for which the input-output relations are well- 
known and we will use a single-resonance Drude- 
Lorentz model for the dielectric permittivity (Sec. ITT|l . 
This, together with some lower bounds on the resonance 
line width, will be sufficient to derive such lower bounds 
that only depend on the frequency of the impinging light 
fSec. EH - We summarize our results in Sec. II VI 



II. THE BEAM SPLITTER MODEL 

We will use a particularly simple model for a beam 
splitter, a single planar slab of thickness I made of a mate- 
rial with a complex refractive index n(cu) = rj{io) + in{u>). 
In this approximation, we consider linearly polarized 
light travelling in a particular direction in space. Fol- 
lowing [l0|, the transmission and reflection coefficients 
T{lo) and R{u), respectively, can be written in the form 
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[n = n(bj)]. Because of absorption, they fulfil |T(w)| 2 + 
\R{uj)\ 2 < 1. 

In the next step, we need to model the complex refrac- 
tive index n(u>). In the linear approximation in which 
the macroscopic polarization responds linearly to an ex- 
ternal electric field, we need to define a complex suscep- 
tibility function x( w )- I n this approximation the quasi- 
excitations associated with the linear response are har- 
monic oscillator-like, hence the susceptibility can only be 
a (possibly infinite) sum of Drude-Lorentz functions, i.e. 
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where lota is the resonance frequencies of the ith res- 
onance with a decay constant or line width 7^, and 
strengths u Pj i that are related to the static dielectric 
permittivity. Note that Drude-Lorentz susceptibilities 
are the only allowed type of response functions. From 
Eq. © it is seen that the susceptibility has a strictly 
positive imaginary part for all positive frequencies. That 
means that even at frequencies far away from all reso- 
nances, a beam splitter will absorb light and we have 
the strict inequality |T(cj)| 2 + |i?(w)| 2 < 1 except for the 
point u> — 0. In particular, for frequencies below all reso- 
nances, we can write for the real and imaginary parts of 
the refractive index to their respective lowest orders in ui 
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This limit is of special importance for our consideration 
since it is the regime in which absorption is lowest while 
n(u>) retains a strong diffractive part. 



III. MINIMAL ABSORPTION PROBABILITY 

Let us assume that we would like to build a beam split- 
ter with a certain splitting ratio x = \T\ 2 /\R\ 2 out of such 
a single slab. Then there are several parameters which 
we can tune: the static permittivity, the frequency of the 
incoming light, the thickness of the slab, and the line 
widths of the resonances. If we restrict ourselves to a 
single-resonance model with a transition frequency wt, 
line width 7, and strength up, we are left with four pa- 
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rameters: the static refractive index 77 = 1 
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the scaled line width 7 = j/ljt, the scaled frequency 
lu = lu/lot, and the scaled thickness d = lutI/c. The 
product rjd is proportional to the optical path length and 
thus determined by the chosen beam splitting ratio x. 
Then, given the frequency u> and the line width 7, we 
minimize the absorption probability p = 1 — |T| 2 — |_R| 2 
over the beam splitter thickness d. For low frequencies 
and narrow line widths we can expand Eqs. Q and (0) 
with respect to those small parameters. The result can 
be cast in the form 
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where a(x) is a numerical coefficient that depends on 
the beam splitting ratio x. FigurcQ]shows the result of a 
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FIG. 1: Coefficient a in Eq. @ for different beam splitting 
ratios x = \T\ 2 /\R\ 2 . The maximum is obtained around a 
symmetric splitting ratio. 

numerical calculation of the coefficient a{x) for different 
x. The maximum of the curve is obtained around the 
symmetric splitting ratio x = 1. At this point, a(l) ~ 
0.9. In the extreme cases when either x — > (perfect 
reflection) or x — > 00 (perfect transmission), a{x) tends 
to zero. However, in the case x — > that does not say 
that the absorption probability vanishes. 

The minimal absorption probability JfJJ still depends 
on two parameters, the resonance line width 7 and the 
frequency of the incoming light, Q. In order to find the 
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absolute minimum of Eq. @ we need to give a lower 
bound on the scaled line width 7. For a non- vanishing 
susceptibility we need non-zero dipole moment which im- 
plies the existence of some (electronic) excited states that 
can be populated by optical excitation. This is because 
the dipole moment is proportional to the off-diagonal ma- 
trix element of the position operator. These electronic 
states are coupled to the electromagnetic vacuum and are 
thus subject to spontaneous decay. Disregarding fabrica- 
tion errors and impurities, spontaneous decay is the lim- 
iting factor when determining the line width. This funda- 
mental process is also responsible for giving lower bounds 
on factorization times in quantum computers 

Because each individual dipole (atom) is surrounded 
by many other dipolcs of the same kind, it will feel not 
only the vacuum fluctuations but enhanced fluctuations 
due to the presence of the surroundings. These effects are 
known as local-field corrections, and several models exist 
to calculate them. In particular, the so-called real-cavity 
model has attracted much attention because it consis- 
tently places the radiating dipoles in vacuum 0, ^ 
rather than in front of a smeared dielectric background 
which is known to cause spurious divergences |l4l flEl Il6| . 

Because we are interested in the low-frequency limit 
in which absorption is very small, we can ignore its in- 
fluence on spontaneous decay and approximate the effect 
of the surroundingatoms by the (classical) real-cavity 
correction factor [l2l Il3| , 
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zero-frequency susceptibility can be written as 
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In that way, we can relate the dipole moment to the static 
refractive index rj as 



d 
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Inserting Eq. (|llfl into JSJl and subsequently into Eq. J7J), 
we find that 



(12) 



where Vt is the volume of a cube with side length equal 
to the transition wavelength At- The product riVr is 
therefore the number of atoms within this cube. 

Finally, we associate Eq. (|12fl with the lower bound 
on the line width 7 of the response function. Hence, we 
obtain the following formula for the minimal absorption 
probability as 



4-7r 2 a(x) 
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Note that the static refractive index 77 appearing in 
Eq. (|13|) has to obtained from the minimization proce- 
dure over the beam splitter thickness d that led to the 
numerical coefficient a(x) in Fig. ^ The dependence of 
the static permittivity e s = 1 + x(0) = rj 2 on the beam 
splitter ratio x is shown in Fig. |21 



where Tq is the free-space decay rate which is defined as 
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where d is the (matrix element of the) dipole moment of 
the transition. The dipole moment is related to the static 
refractive index in the following way. Consider a two- 
level atom with a ground state \g) and an excited state |e) 
(there is no need to consider more complicated situations 
because the single-resonance model corresponds directly 
to the response of an ensemble of two- level atoms). The 
atomic polarizability in the zero-frequency limit can be 
written as |17| 
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where the population numbers of ground-state 

and excited-state atoms, respectively, and n the to- 
tal number density of atoms in the volume V. The 
matrix elements of the dipole operator are denoted as 
di = (e\di\g). The approximation made on the rhs of 
Eq. JHJ is justified because at room temperature and op- 
tical frequencies the thermal occupation of the excited 
state is negligble, and in the linear-response approxima- 
tion there are no population inversions. Therefore, the 
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FIG. 2: Static permittivity e s leading to the minimal value 
Pmin [Eq. @] for different beam splitter ratios x — jT| 2 /|7?| 2 . 

As an example, let us consider the symmetric beam 
splitter with x = 1. The numerical minimization leads 
to a w 0.9 at e s w 6.2. Let us further assume that the 
resonance frequency lot lies somewhere in the UV such 
that we can take nVr ~ 10 9 . With these numbers, the 
minimal absorption probability reads 



10" 6 w 4 . 



(14) 



If additionally the incident light has a frequency that 
lies in the IR region, i.e. Co w 0.1, then we obtain 
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Pmin ~ 10~ 10 . Currently, the best beam splitters that 
have been fabricated show an absorption probability of 
around 2 ppm 0] . This means that on the one hand that 
there is still room for improvement by eliminating techni- 
cal imperfections such as impurities or diffuse scattering 
from interfaces. On the other hand, the lower bound is 
large enough to become an obstacle in future generations 
of beam splitter fabrication. 

Equation l|13|> shows that the minimal absorption prob- 
ability strongly depends on the static permittivity (or, 
equivalently, the static refractive index). For large val- 

r ■ 3/2 3 

ues of rj, Pmin increases as p m in oc e s or p m in oc rj , re- 
spectively. From Fig.[5]it follows that this is the case for 
small beam splitting ratios, i.e. when the slab is highly 
reflecting. For example, choosing x = 0.05 and all other 
parameters as before, we obtain p m ; n « 2 • 10~ 5 ci> 4 , an 
increase by a factor of 20. 

IV. CONCLUSIONS 

We have shown that one can obtain lower bounds on 
the absorption probability of a beam splitter by causal- 
ity arguments. This shows the intimate relation between 
manipulability of light — described by the real part of 
the linear response function — and decoherence or ab- 
sorption which is described by the imaginary part of the 
response function. Although it is intuitively clear that 
such relations must exist, we have presented in a par- 
ticularly simple example how they emerge in a realistic 
physical situation. 

From the Kramers-Kronig relations one infers that a 
beam splitter material which necessarily has r)(uj) ^= 1 
shows absorption somewhere in the frequency spectrum. 
However, it is only the assumption of linear response 



and the resulting validity of Drude-Lorentz models to 
describe the dielectric susceptibility that leads to the ab- 
sorption coefficient being strictly positive for positive fre- 
quencies. In the single-resonance Drude-Lorentz model 
employed here the width of the resonance is limited by 
the spontaneous decay rate between excited and ground 
state. This is certainly a lower limit when all line- 
broadening effects are disregarded. 

The final result, Eq. (| L 3 p . is seen to depend only on the 
transverse resonance frequency lot and the frequency of 
the incident light to. The other parameters (beam split- 
ter thickness static susceptibility e s , and line width 7) 
are fixed by the beam splitter condition (Z), the mini- 
mization (e s ), and the spontaneous decay rate (7), re- 
spectively. For reasonable assumptions on the trans- 
verse resonance frequency, we obtain that the absorp- 
tion limit for a symmetric beam splitter is approximately 
Pmin ~ 10~ 6 tl> 4 . This result will set a possible limit on 
experimental performances where ultra-low absorption 
is necessary, such as in gravitational-wave interferome- 
try. Currently, absorption probabilities of 2 ppm have 
already been achieved, and with progress being made in 
eliminating material impurities, surface roughnesses, and 
other experimental imperfections it is reasonable to as- 
sume that the theoretical limits of the kind that have 
been considered here will impose ultimate performance 
limitations that have to be taken seriously. 
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